Exponentially rapid decoherence of quantum chaotic systems 
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We use a recent result to show that the rate of loss of coherence of a quantum system increases 
with increasing system phase space structure and that a chaotic quantal system in the semiclassical 
limit decoheres exponentially with rate 2A2, where A2 is a generalized Lyapunov exponent. As a 
result, for example, the dephasing time for classically chaotic systems goes to infinity logarithmically 
with the temperature, in accord with recent experimental results. 

PACS numbers: 05.45.+b,03.65.Sq,03.65.Bz,05.40+j 

Decoherence jj],2 which is, in general, the loss of interference effects in a quantal system due to its interaction with 
t-H ■ the environment [3|, is important in fields that rely on the maintenance of delicate quantal coherence, including the 
rapidly-growing areas of coherent control Q], quantum computation || and charge transport in mesoscopic devices 
It also plays a central role in treatments of classical-quantum correspondence |lL2|. Preliminary studies have 
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suggested that systems which are chaotic in the classical limit decohere rapidly @-|]. However, the rate of this 
decoherence and its dependence on system properties is unclear. Here we derive the functional form for the rate 
of decoherence, provide an explicit link between the rate and the extent of the system's phase space structure and 
demonstrate that decoherence occurs exponentially, with a rate determined by a generalized Lyapunov exponent, for 
quantum systems with a chaotic classical analogue. In addition, we note that the decoherence rate is reduced for h 
• sufficiently large, although in general the decoherence of a quantal system may be maximally rapid for a finite value of 
t-H \ h. Finally, we show that our results explain the recent experimental result || that the decoherence rate in mesoscopic 
devices saturates as the temperature T — > 0. All of these results follow from an analysis of decoherence in light of a 
recent approach |U)[ which connects the exponentially fast growth of structure in the distribution dynamics of chaotic 
■ systems to a generalized Lyapunov exponent A2 of the classical dynamics in phase space. 

\ Consider the Wigner representation (= p w ) jllj of a quantal density matrix; it is suitably normalized (Tr[p w ] = 1), 
but is in general not a 'pure state' (Tr^p^) 2 ] 7^ 1). The relevant decoherence equation for p w , derived from the 
quantum system evolving under the potential V(q) and coupled to an external environment through the position 
operator, is IJUl 

S : Q P W m w x , xr ^"t- 1 )" d2n+1 ,rr, d 2 " +1 w n d 2 w 

jS; ~W = {H ^ > + E 2 2n (2n + 1) ,a g2w+ i ng) g^P +D — p . (1) 

\ The first term on the right of this equation is the Poisson bracket or classical Liouville operator, which generates 
the classical evolution for p w ; the terms in h add the quantal evolution, with the expansion being formally valid 
for a potential V(q) analytic in q. The decoherence is contained in the term dependent on D; its precise form and 
parameter-dependence rests upon assumptions about the form of the system-environment coupling and the spectrum 
of the environment. The standard form Q, for example, has D = 2^MHbT where M is the system mass, T the 
temperature of the environment, ks is the Boltzmann constant and 7 an unknown measure of the system-environment 
coupling. Different forms for D may also obtain |§,|l2|, although D always depends monotonically on the coupling 
and the temperature of the environment. The particular form of the coupling, D-^p w comes from the choice of 
coupling through position; coupling through momentum would give DJ^p w . For full generality it is appropriate to 

use DV 2 p w where V 2 is the Laplacian operator in phase space. By using the relationship ^ = ^ + ^-[A,H] for 
the time-evolution of an arbitrary operator A, Eq. (Q) can also be written as 

We use the V 2 form hereafter; the arguments are entirely unchanged if either the momentum or the position coupling 
is used. 

Our criterion for monitoring the degree of system decoherence is the Renyi entropy @Jl^l S = ln(Tr[(p w/ ) 2 ]). This 
function maximizes at zero for a pure state and is a direct measure of the degree of mixing of the quantum state 
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JHfLU . The time-dependence of S thus provides information about the rate at which a "coherent" quantum pure state 
is transformed into an "incoherent" or statistical mixture of states. It is straightforward to show, using Eq. (|^) that 

The quantity \ p0| (which is in general time-dependent) affords insight insofar as it is the root-mean-square Fourier 
radius of the distribution, and a measure of phase space structure. That is, if we Fourier expand any distribution p 
as p(p, q) — J J dp, dv e in% (w+ v i) p(u, v) where the tilde now represents the Fourier-transformed function, we get that 

2 = Tr[pV 2 p] = Tr[ |Vp| 2 ] = ^ J dk fc 2 |p(fc)| 2 
X Tr[p 2 ] Tr[p 2 ] J dk |p(fc)| 2 1 1 

where the second equality follows from integration by parts and where k = {(J,,f). Since higher Fourier modes are 
related to increased structure, x measures the structure in the distribution, and Eq. (^|) shows that the decoherence 
rate of a quantal system is directly proportional to the extent of system structure in phase space [ [i7)[ . Note that since 
^ is negative-definite, the system can not "recohere" jL6j and hence much of the decoherence is determined by the 
short-time behavior of x(t)- Therefore, as long as the term in D is not initially dominant, we can make quantitative 
predictions of dS/dt based upon the time-dependence of the phase space structure from the Hamiltonian part of the 
evolution for p w Jl7[ . 

Consider first time-independent states (e.g., a mixture of Hamiltonian eigenstates). To first order, i.e. before 
decoherence becomes significant, we can estimate x 2 from the initial p w , and find that lower energy states, with less 
structure, decohere more slowly. Thus, for example, for the nth eigenstate of the harmonic oscillator % 2 oc n so that 
decoherence is proportional to n, explaining the numerical results of Knight and Garraway |jl8|| . Hence we deduce 
that in general stationary distributions with lower energy are more stable with respect to environmental perturbation. 

For time-evolving distributions we need to consider Eq. (||). For nonchaotic systems some insight results from a 
simple analysis since in this instance an adequate measure of the time-dependent x 2 is the time-averaged x 2 ■ For pure 
states of the harmonic oscillator for example, working within the approximation that an initially pure state remains 
so, this criterion yields that minimum uncertainty coherent state are the most stable |Q. 

For chaotic systems insight obtains from examining the behavior of x 2 m the semiclassical limit. In this limit p w is 
known jl0| to behave like the classical result, at least for early times. Consider then the behavior of x 2 for a classically 
chaotic system. Let the equations of motion of a point in phase space be x = f(x). The equations of motion for the 
vectors in the tangent space are obtained by linearizing the dynamics x — ► x + ? around a fiducial trajectory to yield 

where the Jacobian matrix M. has the elements My = at the point x(t). It is straightforward to show |l(| that 

±V P =-MVp, (6) 

that is, the equations for evolution of c and Vp are identical except for a minus sign, where Vp is the phase-space 
gradient of the classical probability density p. This implies that the usual maximal Lyapunov exponent may be 
computed from the distribution gradient evaluated along a trajectory as 

A(x)= lim jH\\7p(x(t))\) (?) 

We now note that x 2 is the distribution average of this gradient. Thus, adaptating the results of the standard 
thermodynamic approach for similar averages which use <r instead of Vp JL3[ , we have p0[ | that 

lim -ln(x) = A 2 . (8) 

t— >oo t 

Here A2, which is independent of trajectories, is a generalized Lyapunov exponent related to the ordinary Lyapunov 
exponent by A2 = A + Q, where Q > measures the local phase-space fluctuations in A (equivalent to a measure of 
the degree of intermittency in chaos |l5|]). The quantity A2 has the particularly simple interpretation M] in that it 
measures the exponentially rapid increase of structure for a chaotically evolving distribution in a Hamiltonian system 
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and is in this sense a fundamental exponent for chaos in distribution dynamics. From this treatment it is clear that 
x(t) grows exponentially rapidly in a chaotic system as 

X 2 W=X 2 (0)[exp(2A 2 t) + C]. (9) 

where C in general fluctuates or grows slower than exponentially with time. 

Thus, substituting Eq. (||) in Eq. ([|) we get that the initial decoherence rate of a quantal system in the semiclassical 
limit, to first order in D is 

^ = -2£x 2 (0)[cxp(2A 2 <) + C], (10) 

an equation which displays a direct correlation between decoherence rates and the generalized Lyapunov exponent of 
the underlying classical system. 

Note that Eq. ( |Io|) implies that in the limit of chaotic semiclassical behavior, all time-dependent quantum states 
decohere exponentially rapidly and hence that there is only a marginal dependence on the initial distribution. This 
is in contrast to non-chaotic systems where, as noted above, the rate is strongly affected by the initial state. 

We have numerically solved the discrete map version of Eq. ([l]), with the full V 2 coupling to the environment, for 
a quantum chaotic Cat Map jjJLEMiJ, to explore the nature of S(t) as a function of a (a scaled U Jn],|2(]]). Results for 
a = 1CP 1 to 10~ 5 are shown in Fig. (1), with D chosen at a typical value of 10~ 6 . Several features are evident: (a) 
in the limit of a — > 0, 5 decreases exponentially rapidly initially with a rate determined by A 2 = A = 0.9624 
consistent with Eq. (|l(i|), (b) the initial decoherence rate increases with a for small a, with maximal falloff for an 
a as 0.001, and (c) for a sufficiently large, the entropy decreases quite slowly. 

Examination of the first quantum correction to Eq. jio] ) provides some insight into the observed dependence of S(t) 
on a. Note first that as h is increased away from and for small t, there is an increase in the amount of structure in 
the distribution. That is, one expects p w to adiabatically follow the classical distribution, but with added fringes p^ ] 
(interference structures). The added structure can be seen, for example, in comparing the contour maps (Figs. (1-3), 
Ref. |2(J) of the evolving distributions. From the first quantal term in Eq. (|l]), we can estimate that the quantum 
correction to an individual Fourier component p(k) is p(k) — > p w (k) as p{k){l + c'h 2 k 3 ) where c' depends on the 
potential and includes all other constants. Thus, for an initial distribution that is not very structured (vanishing 
support at large k), p w is very close to p. However, even when the quantum terms are small, x 2 sums them over all 
the Fourier components, weighted by the factor k 2 . A crude estimate |23|] then yields that the quantum corrected x\ 

relates to the classical estimate Xc as Xq ~ Xc(l + c ^ 2 Xc)- Hence the quantum correction to x 2 can be quite large. 

This analysis of the contribution of an individual Fourier component to the enhancement of structure breaks down 
for K 2 k 3 = 0(1); this happens rapidly, since for classically chaotic dynamics the support for the distribution at larger 
k values increases exponentially fast. Hence, higher-order contributions also have to be estimated, which is in general 
extremely difficult. By this time, however, the second quantal effect enters: Quantum distributions resist the growth 
of structure at scales smaller than h [)22| and the support of the distribution travels slowly across the k as — boundary 
in Fourier space. This is marked by fllOf a slow-down of the growth of x a t C(~)< The interplay between these two 
effects leads to a value of Ti where the quantum system decoheres maximally rapidly for early times. Quantum systems 
for h larger than this critical value are, however, stable against decoherence, as seen in Fig. 1. 

A second interesting effect emerging from this treatment is the dependence of the decoherence rate on D. It has 
been recently observed [0] that the decoherence rate for quantal mesoscopic devices goes to zero slower than any 
power of the temperature T as the systems are cooled to near-zero T . The measure used in these experiments is the 
dephasing time, which is monitored through the change in coherent quantal back-scattering effects in the device (weak 
localization) [@J. In our theoretical analysis, arguing that an increase in the degree of "incoherent" mixing reduces 
this coherent back-scattering effect, a reasonable analogue of the dephasing time may be taken to be the time tu 
at which the system has "sufficiently decohered" , i.e when the entropy has fallen to some fixed fraction of its initial 

value R — 5(0) — 5(i#) = ^( ^[^(fJI)]] ) = 1; f° r example. Since these mesoscopic devices are disordered systems, 
with the electrons experiencing substantial chaotic scattering from impurities, Eq. ( ]To| ) is expected to apply in the 
semiclassical limit. If we now solve for the effect on S(t) of the dominant exponential term from this equation, we get 
the expression 
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(11) 



The 1/D term dominates as T, and hence D,— > and indicates that, independent of the particular power-law 
dependence of D on the temperature, the dephasing time goes to infinity logarithmically with vanishing temperature, 
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i.e, slower than any power-law. Equation ( |ll| ) is, in fact, able to reproduce the general shape of Figs. (1) and (2) of 
Ref. |(|, albeit with free parameters. In addition to this logarithmic dependence on D in the semiclassical limit, a 
different effect shows up as the system becomes more quantal: Since (1) the decoherence rate grows as D\ 2 and (2) 
the growth of structure slows down at \ = 0(h), the decoherence term has a distinctly greater impact if \Hj > 0{fr). 
This behavior can be seen in the slow decrease in S for a = 1CU 2 , 1CU 1 in Fig. 1. 

In summary, we have demonstrated that the decoherence rate of a quantal system coupled to the environment is 
governed by the degree of structure in the system as measured by the quantity x- This has enabled us to make 
various deductions about the stability of distributions under the action of the environment which generalize several 
recent [|jj|7],|8| discussions. We have shown that since the structure of a distribution increases exponentially fast for a 
classically chaotic system, distributions in such systems decohere exponentially fast in the semiclassical limit, with the 
exponent given by the quantity A2. This behavior is manifest, for example, in the observation Q that the dephasing 
time for semiclassical systems which are classical chaotic goes to infinity logarithmically with the temperature. 
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FIG. 1. Time-dependence of the degree of incoherent mixing of quantal distributions as measured by the change in entropy 
(S(t) — S(0)) for the chaotic quantal Cat Map coupled to the environment. The effect of approaching the classical limit is 
shown by varying a from 10 to 10 , where a is a scaled h. D, which is a measure of the coupling to the environment, is 
10~ 6 . See text for details. 
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